AD-A077  356 
UNCLASSIFIED 


MASSACHUSETTS  INST  OF  TECH  LEXINGTON  LINCOLN  LAB  F/G  11/6 

MAGNETIC  MOMENT  VERSUS  TEMPERATURE  CURVES  OF  RARE-EARTH  IRON  GA— ETC (U) 
JUL  79  G  F  DIONNE  F19628-78-C-0002 

TR-534  _  ESD-TR-7Q-179 _ NL 


k. 


nrjrwiaitilfirJ 


ESD-TR-79-179 


Technical  Report 


G.  F.  Dionne 


Magnetic  Moment  Versus 
Temperature  Curves 
of  Rare-Earth  Iron  Garnets 


Prepared  for  the  Department  of  the  Army 
under  Electronic  Systems  Division  Contract  F19628-78-C-0002  by 

Lincoln  Laboratory 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 


Lexington,  Massachusetts 


Approved  for  public  release;  distribution  unlimited, 


Non-Lincoln  Recipients 

PLEASE  DO  NOT  RETURN 

Permission  is  given  to  destroy  this  document 
when  it  is  no  longer  needed. 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
LINCOLN  LABORATORY 


MAGNETIC  MOMENT  VERSUS  TEMPERATURE  CURVES 
OF  RARE-EARTH  IRON  GARNETS 


G.  F.  DIONNE 
Group  33 


TECHNICAL  REPORT  534 

11  JULY  1979 

Approved  for  public  release;  distribution  unlimited. 

MASSACHUSETTS 

79  11  26  103 


LEXINGTON 


ABSTRAC  T 


Molecular -field  coefficients  of  the  heavy  rare-earth  iron  garnets  (Tb,  Dy,  Ho, 
Er,  Tm,  and  YbIG)  have  been  determined  by  fitting  calculated  and  experimental 
magnetic -moment -vs-temperature  curves.  For  these  compositions,  spin  cant¬ 
ing  appears  to  cause  a  significant  reduction  in  the  net  moment  of  the  c  sublattice. 
The  results  are  compared  with  values  of  coefficients  derived  from  paramagnetic 
susceptibility  measurements  and  with  other  values  determined  without  taking  into 
account  c -sublattice  canting.  The  capability  of  the  computer  program  has  been 
expanded  to  permit  magnetization-temperature  calculations  of  compositions  that 
contain  as  many  as  three  different  rare  earth  ions  in  the  c  sublattice.  _ . 
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MAGNETIC  MOMENT  VERSUS  TEMPERATURE  CURVES 
OF  RARE-EARTH  IRON  GARNETS 


I.  INTRODUCTION 

The  magnetic-moment-versus-temperature  (ng-vs-T)  characteristics  of  rare-earth  iron 
garnets  have  been  important  in  microwave  ferrite  development  and  more  recently  in  the  emer¬ 
gence  of  ferrimagnetic  thin-film  technology  for  bubble  memory  devices.  By  determining 
molecular-field  coefficients,  the  curves  of  these  materials  may  be  computed  from  the  Neel 
theory,  and  a  basis  may  be  established  for  predicting  the  curves  of  compositions  with  nonmag¬ 
netic  ions  substituted  for  Fe3+  or  rare-earth  irons. 

In  addition,  it  is  also  possible  to  compute  curves  for  compositions  containing  more  than  one 
type  of  rare-earth  ion.  The  purpose  of  this  report  is  to  present  the  results  oT  investigations 
with  a  computer  program  that  was  used  to  determine  the  molecular-field  coefficients  and  ef¬ 
fective  angular  momentum  quantum  numbers  of  heavy  rare-earth  iron  garnets. 

II.  THEORETICAL  METHODS 

In  work  reported  previously,  the  molecular-field  coefficients  N. .  of  Y-FetO,,  were  deter- 

2  —4  ^3  ^  o  x  Ci 

mined  by  a  trial-and-error  procedure.  To  arrive  at  a  set  of  coefficients,  N. .  were  adjusted 

until  the  computed  ng-vs-T  curves  produced  exact  fits  with  the  experimental  curves  for  the  same 

materials.  For  the  above  two  materials,  it  was  found  that  the  octahedral  (a)  and  tetrahedral  (d) 

intrasublattice  and  intersublattice  coefficients,  Naa,  Ndd,  and  Nad>  originally  determined  for 

Y3Fe5Oi2,  were  equally  suitable  for  Gd^Fe^O^^-  However,  the  effects  of  magnetic  Gd5  +  ions 

in  the  dodecahedral  (c)  sublattice  were  accounted  for  with  values  of  N  ,  and  N  ,  listed  in 

ac  cd 

Table  I.  For  the  purpose  of  establishing  a  good  fit  to  experimental  results,  it  was  convenient 

to  assume  that  the  small  values  of  Ncc  were  negligible  (see  Fig.  1). 

Since  quenching  of  the  orbital  angular  momentum  is  insignificant  in  the  rare-earth  family, 

the  value  of  the  spectroscopic  splitting  factor  g  is  equal  to  the  pure  spin  value  of  2.0  only  in 
3  +  .  .  c 

the  case  of  Gd  (an  S-state  ion).  With  every  other  ion,  it  is  derived  from  the  relation  g  = 

(Lc  +  2Sc)/[Lc  +  Sc]  where  Lc  and  Sc  are  the  orbital  and  spin  quantum  numbers  of  the  ground 
term,  and  should  be  combined  with  the  total  angular  momentum  Jc  and  the  Bohr  magneton  jig 
in  the  Neel  theory.  In  these  cases,  Jc  replaces  Sc  and  the  magnetic  moment  per  ion  becomes 

gc  ^BJc’ 

In  Appendix  A,  the  format  for  the  modified  Neel  theory  is  presented  for  the  general  case 
with  as  many  as  three  different  c-sublattice  rare-earth  ions.  Because  the  intrasublattice  coef¬ 
ficients  Ncc  are  negligible,  it  may  be  assumed  that  any  rare-earth  ion  is  coupled  only  to  the 
iron  in  the  a  and  d  sublattices  and  does  not  interact  any  other  rare-earth  ion.  In  this  way,  it 
is  possible  to  treat  these  cases  as  if  separate  c-sublattices  existed  for  each  type  of  rare-earth 
ion,  with  variables  x,  y,  and  z  designating  the  respective  fractions  of  the  different  rare-earth 
ions  in  the  dodecahedral  sites.  A  complete  listing  of  the  Fortran  program  (compatible  with  the 
IBM/370  computer)  is  given  in  Appendix  B. 

RI.  RESULTS  OF  CALCULATIONS 

With  a  trial-and-error  procedure  identical  to  that  used  in  the  previous  work,  molecular- 
field  coefficients  Nac  and  Ncd  were  determined  for  Tb,  Dy,  Ho,  Er,  Tm,  and  Yb  iron  garnets5'8. 
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Fig.  1.  Comparison  between  theory  and  experiment  for  ng-vs-T  curves 
of  YzGd3_zFe50i2.  using  data  of  Vassiliev  et  al.  (Ref.  13),  Harrison  and 
Hodges  (Ref.  9),  and  Anderson  et  al.  (Ref.  14). 

Examples  of  the  precision  of  the  fits  between  calculated  curves  and  experimental  data  points  are 
given  in  Figs.  2  through  7.  For  these  compositions,  it  was  necessary  to  introduce  an  effective 
angular  momentum  <  Jc  to  produce  an  exact  fit  at  T  =  0  K.  The  values  of  JJ,  required  for  this 
purpose  are  given  in  Table  I,  and  suggest  that  a  canting  effect  or  departure  from  collinearity 
within  the  c  sublattice  likely  accounts  for  the  reduced  c-sublattice  moments.  The  extent  of 
canting  has  been  estimated  by  computing  a  canting  angle  <p,  defined  by  the  relation  JJ,  =  Jc  cos  <p 
(see  Table  I). 

Because  of  the  generally  weak  magnetic  exchange  coupling  associated  with  c  sublattices, 
this  effect  may  be  attributed  to  magnetocrystalline  anisotropy.  This  argument  is  consistent 
with  the  fact  that  the  correction  of  Jc  for  canting  was  not  necessary  with  Gd3+,  the  only  ion  of 
the  group  without  an  orbital  contribution  to  Jc  and  the  spin -orbit  coupling  energy  required  for  a 
first-order  anisotropy  effect. 

To  demonstrate  the  program's  capability  with  different  rare-earth  ions  in  the  same  compo¬ 
sition,  complete  curves  are  presented  in  Fig.  8  for  the  system  Y2>7_xGdxDyQ.3Fe5  °12  together 
with  available  experimental  data/*  and  a  comparison  between  theoretical  and  experimental  compen¬ 
sation  temperatures11'  for  ErxGd3_xFe5Oi2  is  given  in  Fig.  9. 

IV.  DISCUSSION 

In  work  reported  by  Aleonard,11  molecular -field  coefficients  were  determined  for  these 
compositions  from  paramagnetic  susceptibility  measurements  above  the  Curie  temperatures. 
Although  these  values  (see  Table  I)  could  not  be  expected  to  produce  accurate  n^-vs-T  curves, 
it  is  interesting  that  their  relative  magnitudes  agree  at  least  qualitatively  with  those  determined 
in  this  work. 

More  recent  values  of  these  coefficients  have  been  reported  by  Brandle  and  Blank1  ^  and  are 
also  included  in  Table  I  for  comparison.  These  values  were  determined  in  a  manner  similar  to 
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Fig.  2.  Comparison  between  theory  and  experiment 
for  the  ng-vs-T  curve  of  Tb3Fe50t2,  using  data  of 
Pauthenet  (Ref.  5). 


Fig.  3.  Comparison  between  theory  and  experiment 
for  the  ng-vs-T  curve  of  Dy3FesOi2,  using  data  of 
Pauthenet  (Ref.  5)  and  Vassiliev  et  al.  (Ref.  7). 


Fig.  4.  Comparison  between  theory  and  experiment 
for  the  ng-vs-T  curve  of  Ho3FesOi2,  using  data  of 
Pauthenet  (Ref.  5)  and  Vassiliev  et  al.  (Ref.  7). 


Fig.  5.  Comparison  between  theory  and  experiment 
for  the  ne-vs-T  curve  of  Er3FesOi2,  using  data  of 
Pauthenet  (Ref.  5). 
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Fig.  8.  Comparison  between  theory  and  experiment 
for  n3-vs-T  curves  of  Y2.7-xGdxDy0.3FesOi2,  using 
data  of  Harrison  and  Hodges  (Ref.  9). 


Fig.  9.  Comparison  between  theory  and 
experiment  for  compensation  tempera¬ 
tures  of  ErxGd3_xFe5Oi2.  using  data  of 
Villers  et  al.  (Ref.  10). 
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the  procedure  used  here,  but  the  work  did  not,  in  general,  include  a  fit  to  the  low -temperature 
points  at  T  ~  0  K.  Consequently,  the  effect  of  canting  was  not  taken  into  account  and  reduced 
values  of  J  were  not  used  in  the  determination  of  N  and  N  ..  With  these  higher  values  of 

C  aC  LU 

Jc,  the  molecular -field  coefficients  required  to  produced  accurate  n^-vs-T  curves  over  the 
higher  temperature  range  would  necessarily  have  to  be  smaller  than  those  reported  in  this  work. 
The  differences  in  magnitude  of  the  coefficients  are  discernible  in  Table  I  by  inspection  of  the 
respective  values  and  are  consistent  with  this  observation.  However,  it  should  also  be  pointed 
out  that  j  practical  purposes  either  set  of  coefficients  will  produce  accurate  fits  over  the 
higher  temperature  ranges  where  most  applications  are  likely  to  be  found,  provided  that  the 
appropriate  values  of  are  employed  with  each  set. 
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APPENDIX  A 


For  the  three  sublattice  case  with  as  many  as  three  different  magnetic  ions  in  the  c  sublattice, 
the  Neel  theory  must  be  augmented  as  follows: 

M=  |Md-Ma-Mcx-Mcy-Mcz|  ,  (A-t> 

where  the  subscripts  d,  a,  and  c  refer  to  the  tetrahedral,  octahedral,  and  dodecahedral  sub¬ 
lattices,  and  x,  y,  and  z  represent  different  c-sublattice  magnetic  ions.  For  each  sublattice, 
arm  =  ar (0)  Bg  (x.),  with 


Sd^B 

kT 

Saga|1B 

kT 


(N .  .M  ,  +  N,M  +N.M  +N.M  +N.M  ) 
dd  d  ad  a  cdx  cx  cdy  cy  cdz  cz 


(N  ,M ,  +  NM+N  M  +  N  M  +N  M) 
ad  d  aa  a  acx  cx  acy  cy  acz  cz' 


xcx  =  J  CXkT^P  <NcdxMd  +  NacxMa  +  NccxMcx>  • 
xcy  =  J,cygcy^B(NcdyMd  +  NacyMa  +  NccyMcy)  ' 

xcz  =  J’czgcz**B<NcdzMd  +  NaczMd  +  NcczMcz>  *  <A 

where  Bg  is  the  Brillouin  function,  k  is  Boltzmann's  constant,  and  is  the  Bohr  Magneton. 

For  practical  purposes,  it  may  be  assumed  that  N  ,  N  .  and  N  “0. 

cox  c  c  ccz 

At  T  =  0  K,  the  magnetic  moments  per  formula  unit  are  given  by 
Md(°)  =  3gdSdpBN(l  -  kd)  (i  -  0.1  ka)  , 


Ma(°)  =  2gaSa(xBN(1-ka>  d  -  “a  > 


“cx^-^c^'c^B^-V  ' 

Mcy(0)  =  3ygcyJ'cypBN(l-kc)  , 

“cz^-^cz^cz^-V  * 

where  N  is  Avogadro's  number,  and  kd>  ka,  and  kc  are  the  fractions  of  nonmagnetic  ions  in  the 
respective  sublattices,  and  x,  y,  and  z  represent  the  respective  fractions  of  magnetic  ions  in 
the  c  sublattice  (i.e.,  x  +  y  +  z  =  1). 
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o  o  o  o  n  ti 


FILE:  DI3PDB2  FORTRAN  A 


1/08/79  15:07 


B.I.T.  LINCOLN  LABORATOBT 


NAIN  PROGRAM  FOR  3  MOMENTS —  OSING  SUBROUTINE  HOH3  FOB  CON¬ 
VERGENCE.  G.  DIONNE 
IMPLICIT  REAL*B  (A-H, K, N, N, O-Z) 

RE AL*8  JCPX.JCPY, JCPZ 
BE AL*4  TT,  YY 

COMMON  /CHI/  SD, S A, JCPX.JCPY, JCPZ , HUB , K , BN , GD,GA, GCX.GCY ,GCZ, X, I, Z 
COMMON  /CM 2/  NDD,NAA,NAD,NACX,NACY,NACZ,NCDX,NCDY,NCDZ,NCCX,NCCY, 
1NCCZ 

COMMON  /CM  3/  KD, K A, KC 

COMMON  /ANS/  M , MD 1  ,HA 1 , MCX 1 , HCY1 , MCZ1 , NBAN S 

DIMENSION  TEMP (200) ,NB(200) ,KAA(21)  ,KDD (21) ,NBB  (21,200) .HSAVE  (21) , 
1 KCC  (5) 

THE  FOLLOHING  ARRAYS  ARE  USED  IN  PLOTTING  ONLY 
DIMENSION  TITLE (9) ,YY (1700) ,TT(1700) , LC (21) 

Tl  AND  T2  ARE  LABELS  FOR  THE  X  AND  Y  AXES 
DATA  T 1/' TEMP '/,T2/' NB  •/ 

DATA  KCC/. 8D0,.7D0,.6D0,.5D0,.4D0/ 

THE  FOLLOHING  THO  LINES  REFER  TO  ROUTINES  ACCESSIBLE  TO  USERS  OF 
LINCOLN'S  360/67  AND  PERTAIN  ONLY  TO  THE  PLOTTING 
CALL  TEKZZ  (2) 

C  CALL  STOIDVC  G.  DIONNE', 9,0) 

CALL  ENTEBG  (*G.  DIONNE* ,9) 

NGRAPH  =  2 
SD=2. 5 
SA=2. 5 

HUB=9. 27E-21 

K=  1. 38E-1 6 
8N=6.023R*23 
GD=2. 0 
GA=2.0 
HR ITS (6,208) 

208  FORMAT  (•  INPUT  KA,KC,KD') 

READ (5,*)  K A, KC, KD 
HRITE (6,209) 

209  FORMAT  ('  INPUT  X,Y,Z') 

READ  (5,  *)  X ,  Y  ,  Z 
HRITE (6,210) 

210  FORMAT  ('  INPUT  GCX ,GCI , GCZ • ) 

READ  (5 , *)  GCX , GCY , GCZ 
HRITE (6,212) 

212  FORMAT  ( '  INPUT  JCPX.JCPY, .JCPZ') 

R  EAD  (5 ,*)  JCPX.JCPY, JCPZ 
IF  (Y.EQ.O.)  JCPY  =  2. 

IF  (Z. EQ. 0. )  JCPZ  =  2. 

IF (Y.EQ.O.)  GCY  =  2. 

IF  (Z.EQ. 0. )  GCZ  =  2. 

HRITE (6,211) 

211  FORMAT  (•  INPUT  N ACX,N ACY, N ACZ , NCDX, NCDY , NCDZ • ) 

READ (5,*)  N ACX, N AC Y,N ACZ, NCDX,  NCDY, NCDZ 

NCCY  =  0. 

NCCZ  =  0. 

NCCX  -  0. 

C  INITIAL  KA  VALUE 

KA  A  (1)  -.00 

C  INITIAL  KD  VALUE 
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KDD(1)  =.00 

C  I NCREH ENT  FOR  KA  AND  KD 

DK=.01 

C  I NCREMENT  FOR  TEN PEBATORE 

DT=  10. 

C  NUH3ER  OF  CURVES  PER  PAGE  (NUMBER  OF  KD'S  USED) 

ICURVE=1 

DO  59  L= 1 , ICUH VE 
KAA  (L)  =KA A  (1)  f  (L-  1)  ♦  DK 
59  KDD  (L)  =KDD  (1)  ♦  (L-1)  *DK 
DO  218  LRR=1,5 
KC  =  KCC ( L  R  R ) 

CHANGED  1,  6  TO  1,1 
DO  100  LK A= 1,1 
KA  =  KA  A  (LK A) 

THE  FOLLOWING  3  LINES  ARE  RELEVANT  ONLY  TO  THE  PLOTTING 
IXZ  =  0 

NM  AX=— 25. E*25 
NMIN=25.E+25 
DO  200  LK D= 1 , ICURVE 
C  KD=KDD (LKD) 

NDD=-30.4*  ( 1. 87*KA) 

NAA=-65.*(1.-1.26*KD) 

NAD=97.* (1.-. 25*KA-.38*KD) 

TEMP (1) =0. 0 
TEMP (2) =20. 

IF  (KD. LT. .  10.OR.KA.LT.. 10. OR. KC. LT.  .  1  0)  TEMP (2) =4  0. 

IF  (TEMP(2)  .EQ.40.)  GO  TO  21 
TEMP  (29) =293. 

TEMP  (35) =343. 

DO  31  L=2 , 28 

31  TEMP  (L) =TEMP  (2) ♦  (L-2) *DT 
DO  32  L=30 , 34 

32  TEMP(L) =TEMP(28) + (L-29) *DT 
DO  33  L=36 , 72 

33  TEMP  (L) =TEMP{34) +  (L-35) *DT 
GO  TO  22 

21  TEMP(28) =293. 

TEMP (34) =343. 

DO  34  L=2 , 27 

34  TEMP(L)=TEMP(2)  +  (L-2) *DT' 

DO  35  L=29 ,33 

35  TEMP(L) =TEMP(27) ♦ (L-28) *DT 
DO  36  L=35, 72 

36  TEMP(L)=TEMP(33) ♦ (L-34) *DT 

22  CALL  MOM3 
Ml=1. 

M2=-9999. 

DO  70  L= 1 , 7 2 
T=TEMP  (L) 

CALL  CO NT IN (T,ITER,IER) 

IF  (IER.EQ. 1)  GO  TO  150 
IF  (TEH  P ( 2)  . EQ.20..AND.L.EQ. 29)  M1  =  M 
IF  (TEMP (2)  . EQ.20..AND.L.BQ.35)  M2=N 
IF  (TEH P (2) -EQ.40. .AND. L.EQ. 28)  S1  =  H 


r.  n  no 
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IP  (TBH P  (2) .EQ. 40.. AMD. L. EQ. 34)  H2  =  H 
NBB (LK0,L) =NB ANS 
70  CONTINUE 

150  IP  (IER.EQ. 1)  LC  (LKD) *L- 1 
IF(IES.Eg.O)  LC (LKD) =  L 
HSA VE (LKD) -  (H 1-H2)/H 1 
LS-LC (LKD) 

THE  FOLLOWING  DO  LOOP  PORHS  THE  ARRAI  OP  ABSCISSAS  POR  THE  POIHTS 
TO  BE  PLOTTED 
DO  77  LA=1 ,LE 
IXZ=IXZ»1 

IP  (TER P (LA) . EQ. 150.)  1150  =  IXZ 
77  TT (IXZ) =TEHP (LA) 

200  COHTIMOE 

WRITE (6,4)  KC 
4  PORHAT (//'  KC=»,P5.2) 

WRITE (6,888) 

DO  98  LH= 1 , ICORTE 

98  WRITE  (6,999)  KA, KDD (LH) ,NBB (LN, 1)  ,TEHP (LC (LN) ) , HS AVE (LM» 

888  PORHAT (/10X,'KA,,7X,,KD,,7X#,NB  AT  TEHP=0'  ,7X, *TEHP  LI HIT* , 

17X, • (H1-H2)/H1*) 

999  PORRAT (8X,P5. 3, 4X,?S. 3 ,8X,P8. 3, 10X.F8. 3,7X , E12. 5) 

THE  FOLLOWING  DO  LOOP  PORHS  THE  ARRAY  OP  ORDINATES  POR  THE  POINTS 
TO  BE  PLOTTED 
IN=0 

DO  61  LKD=1,IC0RVE 
LP-LC (LKD) 

DO  61  LZ=1,LP 
IN=IN*  1 

IP  (LC  (LKD)  .  Eg.  1 )  GO  TO  61 
IP  (NBB (LKD,LZ) .GT.NHAX)  NHAX=NBB ( LKD, LZ) 

IP  (NBB  (LKD.LZ)  .LT.NHIN)  NHIN=NBB ( LKD, LZ) 

61  YY  (IN) =NBB  (LKD,LZ) 

C  WRITE  OUT  VALUE  OF  NB  AT  150  DEGREES 

WBITE (6, 1001)  YY  (1150) 

1001  PORHAT  ( '  NB  AT  150  DEGREES  =  *,E13.6) 

CALL  TEKBRS ( 1 ) 

CALL  PLOTLN (LE.TT.O. 0,0.0, YY, 0.0, 0.0,0, 'T(DEG.  KELVIN) », 1 4, 
l'NB',2, ’HAGNETIC  HOHENT ' , 1 5, NGRAPH) 

CALL  FILPRG(I) 

100  CONTINUE 
218  CONTINUE 

C  THIS  CALL  SIGNALS  THE  END  OP  PLOTTING 

99  CALL  EXIT 
STOP 

END 

C  CALCULATION  OP  H  VS  TEHP  POR  CARNETS  WITH  GD  (♦03) 

SUBROUTINE  HOH3 

IHPLICIT  REAL*8  { A- H, K, H, N , 0- Z) 

RBAL*8  JCPX,JCPY, JCPZ 

COflHON  /CHI/  SD, SA , JCPX , JCP Y, JCPZ ,HUB,K,BN, GD, G A, GCX,GCY , GCZ , X, Y, Z 
COHHON  /CH 2/  HDD, N A A, N AD, N ACX , NAC Y, NACZ , NC Dt .NCDY ,NCDZ,NCCX ,HCC Y , 
1NCCZ 

COHHON  /CH3/  KD, KA,KC 

COHHON  /ANS/  H,HD1,HA1,HCX1,HCY1,HCZ1,NB 
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PILE:  DI3PUB2  FORTRAN  i  1/08/79  15:07  B.I.T.  LINCOLN  LABORATORY 


P- BOB* BN 

BDO-3. *GD*SD*P* (1.-KD) * ( 1. 1*KA) 

HAO=2. *GA*SA*P* (1 .-KA) *  (1.-KD**(5. «)  > 

HCXO  »  3. 0*X*GCX* JCPX*P* ( 1.-KC) 

HCIO  =  3.0*Y*GCY*JCPT*P*<1.-IEC} 

BCZO  *  3.0*Z*GCZ*JCPZ*P*{1.-KC) 

Cl«1./(2.*SD) 

C2=(2.*SD*1.)*C1 
C3=  1./  (2.  *SA) 

C4=(2.*SA»1.) *C3 
C5X  »  1./(2.*JCPX) 

C5I  *  1./(2.*JCPY) 

C5Z  *  1./(2.*JCPZ) 

C6X  *  (2. *JCPX* 1. ) *C5X 

C6I  *  (2. * JCPY* 1. ) *C5Y 

C6Z  =  (2. *JCPZ* 1. ) *C5Z 

BD=BDO 

HA-BAO 

BCI=BCZO 

BCI=BCYO 

BCZ=BCZO 

BD1-BDO 

BA1-BAO 

BCX1=BCXO 

acii=Bcro 

BCZ1=BCZO 

BETOBN 

C* ********** ************ *************** •**••*••***** **•*•****••* ******** 

ENTBT  CONTIN  (T,ITER,IEB) 

ITBB=0 

IER=0 

IF (T. EQ. 0- 0)  B=DABS (BD-BA-BCX-BCI-BCZ) 

IF  (T.  EQ.  0.  0)  GO  TO  55 
FB  =  BOB/ (K*T) 

IX=1 

9  ITER=ITER* 1 

IF (ITER.GE. 300)  GO  TO  999 

XD=GD*SD*FB* (NDD*BD* N AD*BA*NCDX*BCX* NCDY* BCY*NCDZ*BCZ) 
XA=GA*SA*FR* (N A A*fl A*N AD*BD+N ACX*BCX*N ACY*BCY* B ACZ*BCZ) 
XCX=GCX*JCPX*FR*(BCCX*BCX*NCDX*BD* 

1 HACX*R A) 

XCY=GCY*JCPY*FR* (NCCY*BCY*NCDY*flD* 

1  NAC  Y*B  A) 

C  WRITE  (6,*)  GCZ, JCPZ, FR, NCCZ, BCZ, NCDZ , BD.NACZ.HA 

XCZ=GCZ*JCPZ*FR* (NCCZ*BCZ* 

1 NCDZ*BD*NACZ*BA) 

WRITE  (6,*)  XCZ 
WRITE  (6,*)  XCZ 

BSD=C2*DCOSH (XD*C2) /DSIBH (XD*C2) -Cl*DCOSH (X D*C 1) /DSIBH (XD*C1) 

BS A=C4  *DCOSH ( X A*C4) /DSIMH ( X A*C4) -C3*DCOSH ( XA*C3) /DSINH (XA*C3) 

IP  (X. NE.O. 0)  BSCX=C6X*DCOSH (XCX*C6X) /DSINH  (XCX*C6X) -C5X*DCOSH (XCX* 
1C5X) /DSIBH (XCX*C5X) 

IF  (Y.NE. 0. 0)  BSCY=C6Y*DCOSH  (XCY*C6Y) /DSIBH  (XCY*C6Y) -C5Y*DCOSH (XCY* 
1C5Y) /DSIBH (XC Y*C5 Y) 

IF  (Z- NE.0.0)  3SCZ=C6Z*DCOSH (XCZ*C6Z) /DSINH (XCZ*C6Z) -C5Z*DCOSH (XCZ* 
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FILE:  DI3PUB2  FORTRAN  A 


1/08/79  15:07 


H.I.T.  LINCOLN  LABORATORY 


1C5Z) /DSINH (XCZ*C5Z) 

IF(X.EQ.O.O)  BSCX  =  0.0 
IF(Y.  EQ.  0.0)  BSCY  =  0.0 
IF (Z.EQ.0.0)  BSCZ  =0.0 
BD2=BDO*BSD 
BA2=BAO*BS A 
BCX2=BCXO*BSCX 
BC Y2=BCYO*BSC Y 
BCZ2=BCZO*  BSCZ 
GO  TO  (13,14)  ,IX 

13  IX=2 
BD1=BD 
BA1=BA2 
HCX 1=HCX2 
BCY 1=BCY2 
BCZ 1=BCZ2 
GO  TO  15 

14  IX=1 
BD1=BD2 
H A 1=BA 
BCX1=BCX2 
BCY 1=BCY2 
BCZ1=BCZ2 

15  IF((DABS(BA1-HA) . LT. . 0001. AND. DABS <BD 1-BD) .LT..0001) 

1. AND. (DABS (BCX1-BCX) .LT.. 0001. AND. (DABS (BCY1-BCY) .LT. . 000 1 . AND. 

1  DA BS (BCZ 1-BCZ) .LT.. 0001) )) GOTO  8 
BD=BD1 
BA=BA1 
BCX=BCX1 
BCY=BCY1 
HCZ3RCZ 1 

IP (BD. LT. 0. 0. OR. B A. LT. 0. 0. OR. BCX. LT.O. 0.OR.BCY.LT.0. O.OR.HCZ. LT. 0. 

1 ) GOTO  999 
GO  TO  9 

8  fl=DABS (BD— B A— BCX-BC1-BCZ) 

55  NB=B/(B0B*BH) 

GO  TO  888 
999  IEB=1 
888  RETURN 
END 
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